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Looking for Cosmological Alfven Waves in WMAP Data 

Gang Chen^, Pia Mukherjee^, Tina Kahniashvili^'^, Bharat Ratra^, and Yun Wang^ 

g : ABSTRACT 

Q^l A primordial cosmological magnetic field induces and supports vorticity or 

I Alfven waves, which in turn generate cosmic microwave background (CMB) 

I anisotropies. A homogeneous primordial magnetic field with fixed direction in- 

I duces correlations between the a^.i^m and a;+i_m multipole coefficients of the CMB 

^ I temperature anisotropy field. We discuss the constraints that can be placed on 

[ the strength of such a primordial magnetic field using CMB anisotropy data from 

^ I the WMAP experiment. We place 3 a upper limits on the strength of the mag- 

I netic field of 5 < 15 nC for vector perturbation spectral index n = —5 and 

^ ■ 5 < 1.7 nC for n = -7. 

O '■ 

'-^1 Subject headings: cosmic microwave background — cosmology: observation — 

' ' methods: statistical 

o ■ 
^: 

C/3 . 

^1 1. Introduction 

X ■ 

' The origin of the large-scale part of observed galactic magnetic fields, of ~ yuG (micro- 

■ ' Gauss) strength and apparently coherent over ~ 10 kpc scales, is unknown. They could be 
the consequence of nonlinear amplification of a tiny seed field by galactic dynamo processes. 
An alternate possibility is amplification of a weak seed field through anisotropic protogalac- 
tic collapse and subsequent further amplification via galactic differential rotation. In both 
cases a primordial seed field of strength exceeding 10"^'^ to 10~^^ G, coherent over ~ Mpc 
scales, is apparently needed, and it is often suggested that upto a ~ nG strength seed field 
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might be required. See Kulsrud (1999), Widrow (2002), and Giovannini (2003) for reviews 
of the state of the art in this area. A primordial magnetic field of present strength ~ nG 
can leave observable signatures in the cosmic microwave background (CMB) anisotropy. 

In standard cosmologies vorticity perturbations decay and so do not contribute to CMB 
temperature or polarization anisotropies. The presence of a cosmological magnetic field gen- 
erated during an earlier cpoch^ changes this situation: a primordial magnetic field induces 
and supports vorticity or Alfven waves (Adams et al. 1996; Durrer, Kahniashvili, & Yates 
1998, hereafter DKY). These vector perturbations generate CMB anisotropics.^ The pres- 
ence of a preferred direction due to a homogeneous magnetic field background leads to an m 
dependence of {aim(i*i^), and induces correlations between the ai+i^ra and ai-i^m multipole co- 
efficients of the CMB temperature anisotropy field. Since the CMB anisotropies are observed 
to be random Gaussian^, it is known that such a contribution can only be subdominant. 

We use the observed (ai-i,mai*fi m) correlations measured by WMAP to place constraints 
on the strength of a homogeneous primordial magnetic field. The angular brackets here de- 
note an ensemble average, and the overbar indicates an average over m for each I. Limited by 
cosmic variance uncertainties, this would be the useful measure to characterize the signature 
of a homogeneous primordial magnetic field. 



^Quantum fluctuations during an early epoch of inflation can generate a primordial nG magnetic field, 
coherent over very large scales (see, e.g., Ratra 1992; Bamba & Yokoyama 2004) 

^In the future one can hope to constrain vector modes through their eflcct on CMB anisotropy polarization 
anisotropy spectra. CMB polarization spectra that result from vector perturbations due to a primordial 
magnetic field have been discussed by Seshadri & Subramanian (2001), Pogosian, Vachaspati, & Winitzki 

(2002) , Mack, Kahniashvili, & Kosowsky (2002), and Subramanian, Seshadri, & Barrow (2003), while Lewis 
(2004) considers the case of vector modes supported by free-streaming neutrinos. 

^Colley, Gott, & Park (1996), Mukhcrjec, Hobson, & Lasenby (2000), and Park ct al. (2001) are some 
early discussions of the Gaussianity of the CMB anisotropy. More recent discussions of the Gaussianity of the 
WMAP CMB anisotropy data are in Komatsu et al. (2003), Colley & Gott (2003), Chiang ot al. (2003), Park 
(2004), Eriksen et al. (2004a, 2004b), Coles et al. (2004), Viclva et al. (2004), Copi, Huterer, & Starkman 

(2003) , Hansen et al. (2004), Gurzadyan ct al. (2004), and Mukhcrjec & Wang (2004). The simplest inflation 
models predict Gaussian fluctuations (see, e.g., Fischler, Ratra, & Susskind 1985; Ratra 1985), and this is 
consistent with most observational indications (see, e.g., Peebles & Ratra 2003). While there are indications 
of mild peculiarities in some subsets of the WMAP data, for example the apparent paucity of large-scale 
power (Spcrgcl ct al. 2003; see Gorski et al. 1998a for a similar indication from COBE data) and the 
differences between data from different parts of the sky (see papers cited above), foreground contamination 
(see, e.g.. Park, Park, & Ratra 2002; Mukherjee et al. 2002, 2003; de Oliveira-Costa et al. 2003; Bennett et 
al. 2003b; Tegmark, de Oliveira-Costa, & Hamilton 2003) and other systematics might be responsible for 
part of this. 
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The model on which we base our analysis is introduced in §2. Our analysis of the 
WMAP data and our results are discussed in §3. We conclude in §4. 



2. CMB temperature anisotropies generated by Alfven waves 

We assume that the homogeneous magnetic field B is generated prior to the time of 
recombination. Such a field could be generated during the electroweak phase transition (see, 
e.g., Vachaspati 1991; Sigl, Ohnto, & Jedamzik 1997; Giovaninni & Shaposhnikov 1998), or 
by an a-effect dynamo driven by collective neutrino-plasma interactions (Semikoz & Sokoloff 
2004). The energy density of this field, i?^/(47r), must be small, to prevent a violation of 
the cosmological principle, and so may be treated as a first order perturbation. Accounting 
for the high conductivity of the primordial non-relativistic plasma (with v <^ 1 where v is 
the velocity field of the plasma), we may use the infinite conductivity, frozen-in condition, 
E + V X B = 0. We also assume that charged particles are tightly coupled to the radiation. 
We write B = Bq + Bi, where Bi denotes the first order vector perturbation in the magnetic 
field (where V • Bi = 0), and v = + f2 with Q being the first order vector perturbation in 
the fiuid velocity (where V • f2 = 0) . 

As a consequence of magnetic fiux conservation the field lines in an expanding universe 
are conformally diluted, Bq oc 1/a^, and the Alfven velocity in the photon-baryon plasma 
during the photon dominated epoch when the energy density pr — + pf, — p-y until 
recombination, is va = -Bo/ V 47r(pi? + pr) = 4 x 10~^(So/10~^G), and is time independent. 
Rescaling physical quantities according to the expansion of the universe, the MHD equations 
result in an equation describing Alfven wave propagation with velocity ^^(b • k) = vaP- 
(DKY) 

n = vl{h-kfn. (1) 

Here b = B0/-B0 is the unit vector in the direction of the magnetic field and an ovcrdot 
denotes a derivative with respect to conformal time tj. Choosing only the sine mode, to 
satisfy the initial condition f2(k, 77 = 0) = 0, we have 

n(k,ri) c:tnoVAkpri, |fio| = — |Bi|. (2) 

-Do 

In cq. (1) we have neglected viscosity so the vorticity solution in cq. (2) is applicable only 
on scales bigger than the damping scale. Assuming that the initial vector perturbation is 
generated by a random process, the two-point correlation function of the vorticity field can 
be written as (Pogosian et al. 2002) 

(Q*,(k)Qo,(k')) = kSij - hkj)S{k) + ieijikiA{k)] 5{k - k'). (3) 
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Here eiji is the totally antisymmetric tensor, and the S{k) and A{k) power spectra describe 
the symmetric and helical parts of the two-point correlation function. We assume that the 
spectra S{k) (= |f2o(^)n and A{k) are given by simple power laws of the scale 1/A; on scales 
larger than the perturbation damping scale l/ki), i.e., for k < kfj we take S{k) = Sok'^/k'^'^ 
and A{k) = Aok"^ / k'^'^'^ . Here 5*0 and Aq are dimensionlcss normalization constants with 
5*0 > Aq, and n and m are spectral indexes. The cutoff scale l/kn is the scale below which 
the magnetic field is damped away, due to the non-infinite value of the conductivity. Using 
^dec ~ 0.3 eV and t^ec ~ 10^^ cm, the comoving magnetic field damping wave number at 
decoupling is knit = tdec) ~ 3 x lO^^^cm-^ (DKY). 

The CMB fractional temperature anisotropy, in direction n on the sky, induced by a 
vorticity perturbation (ignoring a possible dipole contribution from the vector perturbation) 
is (DKY) 

AT 

— (n, k) ~ n • fi(k, rj) = n- riow(^^rfec)- (4) 
Decomposing the CMB fractional temperature anisotropy in a spherical harmonic expansion, 

-jri^) ^J2Y1 (^Irnyimin) (5) 

1=2 m=-l 

and using the definition of the power spectrum C/, 

AT AT 1 °° 

(^(n)^(n')) = ^ 5^(2/ + l)GP,(n • n), (6) 



1=2 



where Pi is the Lcgcndrc polynomial, wc obtain, in the isotropic case, Ci = {aj„fiim)- Here 
the angular brackets denote a theoretical (averaging) expectation value over an ensemble of 
statistically identical universes. In Fourier space this expectation value can be replaced by 
integration over all possible wavenumbers, i.e., (...) — > / d^'k/{2'K)^. 

Computing {al^air„i'), it can be shown that the helical part of vorticity does not con- 
tribute (see Pogosian et al. 2002). Hence in what follows we consider only the symmetric 
part of the spectrum. Detailed computation of C;'s for vorticity perturbations are described 
in DKY where it has been shown that the presence of a homogeneous magnetic field induces 
off-diagonal correlations in multipole space, in particular correlations between I and I ± 2 
multipole coefficients. To quantify this we introduce a second power spectrum defined by 

Di{m) = (a;_i_^ai+i,„) = (a;+i,^ai-i,m). (7) 



The two power spectra, Ci{m) and Di{m), depend on the spectral index n, the normalization 
constant Sq, the Alfven velocity va and the perturbation damping wavenumber ko- The 
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power spectra are defined only for n > —7 (the quadrupole diverges at small A; for n < —7), 
and for n > —1 the results are determined by the damping wavenumber ko- The case n = —5 
corresponds to the Harrison- Peebles- Yu-Zel'dovich scale-invariant spectrum [Ci, Di ~ Z^). 
See DKY for a more detailed discussion. 

The non-zero correlation of temperature for unequal Vs has a simple physical explana- 
tion: The presence of a preferred spatial direction, that of the magnetic field Bq, breaks the 
spatial isotropy of the CMB map, leading not only to an m dependence of the correlators, 
but also non-zero off-diagonal (in / space) correlations.^ The temperature perturbation cor- 
relation between two points on the sky depends not only on the angular separation between 
the two points, but also on their orientation with respect to the magnetic field. 

A simple observational test to detect (or constrain) the presence of a homogeneous 
magnetic field in the Universe is based on computing the Di spectrum of CMB anisotropy 
data.^ For this it is useful to introduce the arithmetic mean over m of the two power spectra, 

1 ^ 



Ci 



21 + 



m=—l 



D, 



21 + 



1 ' 

— 5^ {a'i-i,r, 



According to DKY, 

Ci/Di -- 
Ci 



Vdec 



\n 



iknVo) 

r(-^) 
r(-f) 

Vdcc 



-(n+3)^2 



-n 



3T{-n/2y 



-^+^ n < -1 



n < -1. 



n + 1 



n > -1 



(8) 

(9) 
(10) 

(11) 



Using Bi < Bq, we have |OopA;^ < v\ (see eq. [2]). This inequality must hold on all 
scales inside the Hubble radius at decoupling, k > l/rjdec- With the S{k) spectrum definition 
in eq. (3) we therefore get 2So{k / ko)'^'^^ < v\ for l/rjdec < k < ko, implying 



2So{kDVdecP''^'^ 



< 



n < —3, 



(12) 



"^In a recent paper, Bershadskii & Sreenivasan (2004) show that colhsions between Alfven wave packets 
and their cascades could generate arcminute-scale CMB temperature anisotropies, and argue that this is 
consistent with the WMAP data. 

^Other statistics, related to the Di's here, could also provide useful tests (e.g., Hajian & Souradeep 2003). 
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2^0 < v\ n > -3 . (13) 
So for n < —3 the result is independent of the damping wavenumber kn- 

We now estimate an upper Umiting value of /, /c, beyond which our approximation is no 
longer valid. Just like scalar perturbations (Peebles 1980), vector perturbations are affected 
by collisionless damping. Adding a photon drag force term on the right hand side of the 
vorticity equation (1), we can see that there are no oscillations in the vector perturbation 
case, and that damping occurs on scales slightly larger than the damping scale for scalar 
perturbations, when kcrjdec ~ 10 (DKY), corresponding to Ic ~ 500, beyond which our 
approximation breaks down. 

Inserting the limiting values given for Sq in eqs. (12) and (13) in eq. (9), we find 
for n — —5, 

= 9.04 X lO-^^r^ (^^^ , Di^Ci/2.26, (14) 

and, for n — —7, 

= 8.61 X lO-iOr"^ (^^^ , Di = Ci/2.17. (15) 

In this paper, we use n — —5 and n — —7 as two illustrative cases. These two cases are 
interesting as they correspond to a Harrison-Peebles- Yu-Zel'dovich scale-invariant spectrum 
result for the C/'s and A's, and to a possible infiation model primordial vorticity field 
perturbation spectrum, respectively. These two cases also span the range of constraints that 
can be placed on B using this method, in the range —3 > n > —7, i.e., n = —5 gives the 
weakest and n = —7 the strongest constraint on B. With koVo ^ 0.4 x 10^"^, from eqs. (9) 
and (11) B is not meaningfully constrained for n > —3. 



3. Analysis and results 

We use the foreground cleaned Q, V, and, W band co-added WMAP data (Bennett et 
al. 2003a) to determine the ofi'-diagonal correlations. The data are available in the Hcalpix 
format (Gorski, Hivon, & Wandclt 1998b) at resolution A^sidc = 512. For each value of the 
magnetic field amplitude B we generate 5000 simulations of the CMB sky, and each time 
apply the Kp2 Galactic cut mask prior to computing the model D^'s. The expected value 
of the model Di^s, obtained from the mean of these simulations, is then compared with the 
£);'s obtained similarly from the WMAP data, using the statistic. Gonfidence levels on 
the field strength are derived from the resulting likehhood function. 
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Specifically, each simulation is a realization of the CMB with power spectrum Ci given 
by the best fit fiat-A CDM model with power-law primordial power spectrum (Spergel et 
al. 2003), and with Di the same as that predicted by eq. (10) (or more specifically eqs. [14] 
or [15]) for a given value of B. In other words, we generate a/j„'s such that they satisfy 

instead of 

{aij^ai'm') = Sm.,m'Sl,l'Ci. (17) 

The aim^s are generated upto an l^iax of 512 (corresponding to Healpix resolution A^side = 
256; since the WMAP data are expected to contain useful cosmological information upto 
such an /max)- These are then convolved with the beam functions of each of the Q, V, 
and W radiometer channels (8 in all) to produce 8 maps at Healpix resolution Ngide = 512 
(because the noise maps have this resolution). Independent Gaussian noise realizations of 
rms ao/v^A^ from WMAP are added to the maps, and the 8 maps are co-added weighted 
by Aobs/c"o' where the effective number of observations Aobs varies across the sky, and cxo is 
different for each radiometer channel. This is how each simulation is created. Hereafter the 
same analysis procedure that is applied to the data map is applied to each of the simulations. 
This consists of bringing the co-added map down to Healpix resolution Agide — 256 (since 
we mostly use Di^s only upto an Zmax of 300 in the subsequent analysis), and applying the 
Kp2 sky cut prior to computing the D^'s. 

The whole analysis can be repeated for different values of the spectral index n that 
characterizes the spectrum of the magnetic field perturbations. 

Fig. 1 shows the Di's obtained from the WMAP data (crosses) and the median and 68% 
confidence range contours from simulations for two illustrative values of B. The spread of 
about 10~^ mK^ in the 68% confidence contours for l{l + l)Di is consistent with what is 
expected from cosmic variance alone. 

To compare the likehhoods for different B values, we use the diagonal statistic 

300 T)S~\2 

x^-e ' (18) 

1=2 

where is the WMAP data value, and D^i the average and cr; the standard deviation of 
the D^i, both obtained from 5000 model simulations, for each value of B. Note that we do 
not use the full covariance matrix for the DiS but rather just the diagonal terms. This is 
because 5000 simulations are not sufficient to produce a reliably converged full covariance 
matrix for the Dis. As also noted by Eriksen et al. (2004b), for example, even the above 
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test can provide a just comparison between data and simulations. The likelihood is 
proportional to e~^^/^. We calculate the hkelihood for a few different values of B. 

The hkelihood function obtained for the n — —5 case is shown in Fig. 2. After integra- 
tion, we get a 3 (7 confidence upper limit of S < 15 nG. As we can see from Fig. 2, B—0 
G, which corresponds to pure Gaussian primordial fluctuations, is within the 1 a confidence 
range from the peak of the likelihood function (this 1 a range corresponds to a SB of 3.9 
nG). The jaggedness in the likelihood function is from the fluctuation of mean Di^s used 
in calculating (see Fig. 1). With 5000 simulations, the fluctuation of mean D;'s should 
be of the order of a;/V5000. This results in a fluctuation of 0.1 in the values, or a 5% 
uncertainty in the estimated likelihood. This does not much affect the results of our analysis 
(for example for the 3 a hmit on B we look for a Ax^ of 9, which is not very sensitive to a 
5% uncertainty in the estimated likehhoods). 

For the n = —7 case (see Fig. 3), the 3 a confidence limit is i? < 1.7 nG, and again the 
pure Gaussian primordial fluctuation case with 5 = G is not far from the 1 a confldence 
range from the peak. In this case 1 a corresponds to a, 5B — 0.4 nG. More stringent 
limits are obtained in this case as expected (eq. [15] shows that larger D^'s with a stronger 
Z-dependence are predicted by the model for n = —7). 

The conditions of homogeneity and unidirectionality of the primordial magnetic fleld 
may be a better approximation on some scales rather than others. In each of the above cases 
for n, other ranges in I, such as 2 — 100, 101 — 200, 201 — 300, or 2 — 500, did not indicate 
anything qualitatively different, i.e., B — G remains a satisfactory flt. 

4. Conclusions 

We study off-diagonal correlations of the form Di = {ai-i^mO-t+i m) ^^^^ Y^^^ 

WMAP CMB anisotropy data. Such correlations can result from a homogeneous primordial 
magnetic fleld. We do not flnd signiflcant off-diagonal correlations in the data, which appear 
to be satisfactorily flt by a zero primordial magnetic fleld hypothesis. We place 3 a upper 
limits on the strength of the magnetic fleld of S < 15 nG for spectral index n = —5 
and B < 1.7 nG for n = —7. These two cases are interesting as they correspond to a 
Harrison-Peebles- Yu-Zel'dovich scale-invariant spectrum result for the C/'s and -D^'s, and 
to a possible inflation model primordial magnetic fleld perturbation spectrum, respectively. 
These two cases also span the range of constraints that can be placed on B using this 
method. Future GMB anisotropy data should allow for tighter constraints on a primordial 
cosmological magnetic fleld. 
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Fig. 1. — Off-diagonal power spectra obtained from the WMAP data (crosses) and the 
median and 68% confidence contours obtained from model simulations with magnetic field 
strengths B = (solid lines) and 16 (dotted lines) nG, for a magnetic field perturbation 
spectral index n = —5. 
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Fig. 2. — The likelihood as a function of the strength of the magnetic field for a magnetic 
field perturbation spectral index n — —5. 
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Fig. 3. — The likelihood as a function of the strength of the magnetic field for a magnetic 
field perturbation spectral index n — —7. 



